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Abstract 

Models of spin glasses are studied with a phase transition discontinuous in the 
Parisi order parameter. It is assumed that the leading order corrections to the ther- 
modynamic limit of the high temperature free energy are due to the existence of a 
metastable saddle point in the replica formalism. An ansatz is made on the form of 
the metastable point and its contribution to the free energy is calculated. The Ran- 
dom Energy Model is considered along with the p-spin and the p-state Potts Models 
in their p < oo expansion. 
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1 Introduction 



The mean field replica theory has revealed to be a very powerful method to study spin 
glass models The main problem with the replica method is that it necessarily requires 
an ansatz on the solution. The Random Energy Model (R.E.M.) is a very simplified spin 
glass model which can also be solved without the utilisation of the replicas. It is, therefore, 
a very good testing ground for the replica method since, although very simple, it presents 
the typical features of spin glasses such as replica symmetry breaking and non-ergodicity. 
For T greater than a critical temperature Tc the R.E.M. presents no breaking of the 
permutation symmetry among the replicas which, in absence of external fields, tend to 
have the lowest possible value of the mean overlap. When T <Tc the system undergoes a 
phase transition towards a one-replica-symmetry-broken phase. The corresponding latent 
heat is zero and so the transition is second order. The R.E.M. was first introduced and 
solved by B. Derrida who was also able to calculate the finite-size corrections to the high 
temperature free energy Q . 

In the present work this last result will be obtained by making use of the replica 
method. In doing so, a metastable saddle point in the replica space will be individuated 
as the responsible for the leading order finite-size corrections. The p-spin model and the 
p-state Potts model will also be considered. Both these models tend to a R.E.M. in the 
limit p — > oo. 

The scheme of this paper is as follows: Section 1 is based on reference the Random 
Energy Model and Derrida's results for the finite-size corrections will be presented. In 
section 2 the same results as the contribution of a metastable saddle point will be recovered. 
Section 3 will introduce the p-spin model and repeat the calculation performed for the 
R.E.M. in the formalism of a (p = oo)-spin model. In section 4 the results will be extended 
to the case p < oo. In section 5 the p-state Potts model will be introduced and the p = oo 
limit will be studied. Finally, a. p < oo expansion for the Potts model will be formulated 
and the calculations will be extended to this case. 



2 The Random Energy Model 

It is worthwhile recalling the main results on the R.E.M. in order to establish the notations. 
The model describes the behaviour of any system with a fixed number of energy levels 
whose energies are independently distributed according to a gaussian law. If 2^ is the 
number of levels, the model is defined by the properties: 
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P{E,,Ej) = P{E,)P{Ej). (2) 

The solution can be easily derived using a microcanonical argument. For the free energy 
one finds: 
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where Eq = NJy/\n2 and Pc = 2\/ln2/ J. 

The R.E.M. can also be solved making use of the replicas. If one computes the n-th 
power of the partition function one obtains 
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After averaging over the disorder one has 
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where Vp is the number of pi that are equal to p. The Vp verify the conditions 
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To obtain equation (0) it has been used the fact that for large enough it can be written 
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One can find that for T > y/nTc the dominant contribution of (^) is obtained by taking 



vi = n 
yp>2 = 0. 
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The correspondent expression for Z" is linear in n, so one can use the well known formula 

2^x1 _ 2 
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n-+0 n 

to calculate the high temperature free energy (|^(a)). The low temperature expression of 
@ is instead obtained by taking v^-^p = and Up=p = n/p, where p comes out to be equal 
to T/Te. 

Finally one can write, without deriving it, the expression for the finite size high 
temperature free energy. 
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V2A: - 1 Tc <T < VWTTTc. 
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Derrida derives expression (|12|) for k = 1,2, 3. He also makes the more general guess that 
it should be true for every k > 1. 



3 The Metastable Point 

In this section the replica method shall be used to try to derive equation (12). It has been 
asserted that, for T > Tc, the dominant contribution to the equation (^) is given by the 
choice (p!o|). It will now be considered the effect of only one grouping of m replicas. This 
can be done by taking 



n — m, 
1. 



For large N equation (0) can thus be written as follows 
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where dom stands for 'dominant' and sub stands for 'subdominant'. Z^'dom is given by 
the choice (10) in equation (^) while Z^sub is given by the choice (|l^). One then has 
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where in Z^-^^b all the integer m greater than m = 1 are summed over. One then has 
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The sum in ( p^q ) can be written as an integral in the complex plane over the circuit (C) as 
shown in figure 
One has 
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Both the sum and the integral are not well defined. They can be defined by deforming 
the circuit C into a vertical path as indicated in figure 2. It can then be written 
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The integral on the right hand side of the previous equation can be solved with the aid of 
the saddle point method. Two remarks have to be made on this subject. The first is that 
since the integration is over imaginary dm, and, since the exponent is quadratic in m, the 
saddle point m^p is given by the minimum and not by the maximum in m. Therefore 
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The second remark is that to the saddle point contribution one has to add the residues of 
the integrand that fall on the left of the nisp if the phase constant circuit (|) shifts on its 
right hand side when one makes it pass trough m^p. 

The evaluation of the integrand on m^p gives exactly the last term in equation (|l2|). The 
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are reproduced by the residues that are added. The number k of 
those terms is given by the condition 2 < k < rngp which is equivalent to the condition 
^J2k -ITc < T < VWTTTc in equation (P). The ansatz (|l3|) seems therefore to be 
true, at least in the range of temperatures Tc < T < \/7Tc, where expression ( p^ has 
been proved correct by Derrida. Thus, the choice (|l^) (0) reproduces exactly expression 
(p!^) while other possible choices of the Up seem to lead to lower order contributions. In 
reference |§ this assertion could be proved rigorously for choices of the Up such as 
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= 1 with m > 0, 

= 1 with r > m. (20) 

It can also be argued that this is hkely to happen in general for more complicated choices 
of the Vp where the insertion of more groupings is allowed. Furthermore, these arguments 
do not seem to be dependent on k. If this is so, the choice ( p!3[ ) provides the leading 
order corrections to the high temperature free energy for all T > Tc and expression ([l^ ) 
is correct even for A; > 3. 



4 The j)-Spin Model 

All this may now be extended to the p-spin model which is defined by the Hamiltonian 
[Anp{{s}) = - Jii,j2,...,ij,si---Sp + /i^Si, (21) 

{l<ii<i2<---<ip<N) i 

where h is an external magnetic field and the Jii,i2,---,ip are random variables that obey 
the gaussian law 
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For sake of simplicity, from now on it will be assumed J = 1. If one indicates with Ei the 
energy relative to configuration {cr}i, it is a well known result that 
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where qij indicates the overlap between {{cr}i) and ({o"}j)- In the limit p — > oo then, the 
p-spin model reduces to R.E.M.. 
The expression for Z" is 



= ^ exp 
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where Q,,{s) ^ j^^f s'^sl 



Above Tc = l/2vln2 the model presents no breaking of the replica symmetry. After 
the elimination of all the requested Lagrange multipliers, the high temperature solution is 
a n X n matrix of the form 
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Qab = qo, for a^b, 

Qaa — 0, 



(25) 



For T <Tc the matrix Qab has two parameters qo and qi . In the hmit p ^ oo the solution 
is go = and qi = I and one recovers expression (|3|). For simphcity it has been put h = 
because it will not affect the point. 

Assuming now T > Tc, one can proceed equivalently to what has been done in the 
previous section and evaluate Z"- on a one block matrix Qab of the form 

qi qi qo qo qo qo \ 

qi qi qo qo qo qo 

qi qo qo qo qo 

qi qo qo qo qo /2g) 

qo qo qo qo qo 

qo qo qo o qo qo 
qo qo qo qo o qo 
qo qo qo qo qo o / 

QiBiock is a n X n matrix with an m x m block of elements qi. To represent QiBiock it has 
been chosen n = 8 and m = 4. It is easy to see that 
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which is the same result found in the previous section. Ansatz (13) is therefore equivalent 
to a matrix of the form (p^). It is worth remarking that to obtain Z^iBiock one has to 
sum over all the possible ways of inserting an m x m block in an n x n matrix. 

This result will now be extended to the case p < oo in the fashion of |6|. Before that, 
it is appropriate to briefly show the low temperature behaviour of the model. Under Tc, 
the p < oo expansion leads to the mean field equations 



qo = 0, 



where it has been set 
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One also has 
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Equation (^) for the critical temperature comes from the condition that the value of the 
break point in the order parameter function that maximises the free energy is 1. This is 
because of the nature of the transition which, though second order in the thermodynamic 
sense, is discontinuous from the order parameter point of view [^. 

Assuming now to set T > Tc, one can calculate Z'^suh with the insertion of a block 
as it was done in the p ^ oo limit. Defining T^°° = l/(2\/ln2), the finite-p equivalent of 
expression ( [l2|) can be easily obtained 
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where 
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One also has 
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It is worthwhile noting that, assuming rrigp = 1 in equation ( |35| ) one obtains the equation 
(pO|) for the critical temperature. This means that one can find Tc as the temperature at 
which the block inserted in the metastable matrix disappears and QiBiock coincides with 
the stable saddle point matrix of Z^dom- Furthermore, equation (|28|) for qi is recovered 
as a saddle point equation for Z^sub- In absence of magnetic field, it has been assumed 
'70 = 0. 
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5 The Potts Model 



The results obtained up to now will be extended in this section to the Potts model which 
is defined by the Hamiltonian 
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where p{i) = 0, 1, . . . ,p — 1, and the Jij are as usual random variables obeying a gaussian 
distribution with variance and mean value Jq/N. 

For every p > 4, the model undergoes a phase transition of the same kind as the one 
observed in R.E.M., from a replica symmetric phase to a one-symmetry-broken phase. In 
this work all the ferromagnetic order parameters will be neglected and the focus will be 
put only on the glassy aspects of the model. Hence, the free energy is 
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According to the conventional wisdom, one can assume that in the high temperature phase 
one has go = 0. Therefore the free energy is 
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and the entropy 
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which becomes negative \i T < = \J ^jj^ . Under Tc the solution is given by a one- 
symmetry-broken matrix of elements q°^^ . 

To calculate the low temperature free energy it is useful to use the p vectors e°- (a = 
1, • • • ,|)), defined by the relations 



Cj^e- = p5a,fe - 1 i = l,---,p-l, 



(42) 



where repeated indexes are summed. If = e gi = g, with a little algebra, one gets to 
the low temperature free energy expression 
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where z = (zi, Z2, • ■ ■ ; -^p-i) and k = 

In the hmit p — > oo, the last integral can be exactly solved 
found to be given by the mean field equations 



and the solution is 
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This solution describes a states Random Energy Model with variance proportional to 
{p — 1). The finite size corrections for this limit are already known, so the more general 
case in which the finite-p corrections are included shall be treated directly. 
One can define 
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Recent work has been done on the numerical estimation of the p-dimensional integral 
^(m) Q. Correcting the integral in the free energy expression one obtains the mean field 
equations 
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Setting g = 1 — e, one finds an equation for the critical temperature 
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The finite size corrections are obtained proceeding in the same way as for the p-spin. 
Defining 



one gets 
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and finally 
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In principle, all this is equivalent to what has been done for the ^?-spin and all the consid- 
erations made at the end of section 4 can be repeated here. 



6 Conclusion 

For the high temperature phase of the Random Energy Model, the partition function was 
evaluated on a metastable point that was introduced in order to account for the probability 
that a group of m replicas freezes in a phase which resembles the low temperature one. 
In this way the finite size corrections to the free energy were calculated. The result 
was checked with the one obtained by B. Derrida without the use of replicas. The two 
approaches are totally independent. Derrida proves his result to be true only in the range of 
temperatures Tg < T < V^T^. The results obtained in this work coincide with Derrida's, 
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in particular they coincide for Tc < T < \flTc. The rehabihty of this method does not 
seem to depend on the temperature, provided that T > Tc. Therefore, the equivalence of 
the two results in the range where formula ( [T2| ) can be proved to be true, seems to indicate 
the reliability of the result for all temperatures. In extending this ansatz to the p-spin and 
Potts models it was possible to identify a one-block matrix as the metastable point. Mean 
field equations give, for the elements of the block, the same value as the low temperature 
mean overlap. Furthermore, a p < oo expansion was performed for these models in order 
to extend the results to finite-p p-spin and Potts models. 
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